Hydro dynamically enforced entropic trapping of Brownian particles 
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We study the transport of Brownian particles through a corrugated channel caused by a force 
field containing curl-free (scalar potential) and divergence-free (vector potential) parts. We develop 
a generalized Fick-Jacobs approach leading to an effective one-dimensional description involving the 
potential of mean force. As an application, the interplay of a pressure-driven flow and an oppositely 
oriented constant bias is considered. We show that for certain parameters, the particle diffusion is 
significantly suppressed via the property of hyrodynamically enforced entropic particle trapping. 
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Effective control of mass and charge transport at mi- 
croscale level is in the limelight of widespread timely ac- 
tivities in different contexts. Such endeavors involve Lab- 
on-chip techniques [U12], molecular sieves [3H5], biolog- 
ical [H] and designed nanopores [7], chromatography or, 
more generally, separation techniques of size-dispersed 
particles on micro- or even nanoscales [HHTU] , to name but 
a few. Particle separation techniques use the fact that the 
response of the particles to external stimulus, such as gra- 
dients or fields, depends on their physical properties like 
surface charges, magnetization, size or shape. The neces- 
sary force acting on a suspended particle can be exerted, 
for example, by surrounding walls by neighboring 
particles and molecules via hydrodynamic interactions 
[T2*HT3] . or by external electric fields causing electroos- 
motic flows [15] , electrophoresis [16l [17] , induced-charge 
electrokinetic flows [18H20] , magneto- and dielectrophorc- 
sis [2T] [52], etc. The progress in experiments has trig- 
gered theoretical activities and led to the development 
of the Fick-Jacobs (FJ) approach [53"H2T] , in which the 
elimination of equilibrated transverse degrees of freedom 
provides an effective description for diffusive transport 
along the longitudinal coordinate only. Thus far, this 
approach has mainly been limited to energetic potentials 
generating conservative forces on the particles, as given 
by the first term in Eq. below. 

With this Letter, we overcome this restriction by ex- 
tending the FJ-formalism to the most general force field 
F(r) exerted on particles, which can be decomposed into 
a curl- free part (scalar potential 3>(r)) and a divergence- 
free part (vector potential *(r)), which constitute the 
two components of the Helmholtz's decomposition theo- 
rem, 



F(r) = -V$(r) + V x *(r) 



(1) 



As an application that admits a simple interpretation of 
the divergence-free force, we consider particle transport 
caused by the interplay of a pressure-driven flow [2Blf3"3"] 
and a constant bias acting in longitudinal channel direc- 
tion (here x) [34lf37] , This in turn then yields our major 
finding; namely the phenomenon of hydrodynamically en- 




p(r)+Ap 



Figure 1. (Color online) A segment of the corrugated channel 
at cross section z = 0, confining the overdamped motion of 
point-like Brownian particles. The wiggling profiles are given 
by periodic functions to±(x), see Eq. (Ill; a unit cell is marked 
by the separating dashed lines. The quantities Acj and AQ 
denote the minimal and maximal channel widths, respectively. 
The (blue) field lines depict a typical external force field F(r), 
acting on the particles which are the result of a competition 
of a constant bias / and an oppositely oriented divergence- 
free force from the driven solvent induced by a pressure p(r) 
with pressure change Ap along a unit cell. This results both 
in vortices and stagnation points, note the solid (red) circles. 



forced entropic trapping (HEET) of Brownian particles. 

We start by considering spherical Brownian particles of 
radius R suspended in a solvent of density p and dynamic 
viscosity rj. The latter fills a planar, three-dimensional 
channel with confining periodic walls at y — w±(x), with 
period L, and plane walls placed at z — and z = H, 
with H L, see in Fig. [lj Assuming that (i) the particle 
suspension is dilute and (ii) the particles are small, i.e., 
R <C A.lo with the density comparable with p, implies 
that inertial effects, hydrodynamic particle-particle and 
particle-wall interactions, and, as well, effects initiated 
by rotation of particles can safely be neglected [12] [38] . 

Passing to re-scaled variables r — > rL, $ — > $k B T, 
* -> * k B T, F -> Fk B T/L, and t->rt, in terms of the 
relaxation time r = 6iri]RL 2 / (k B T), where r = (x, y, z) T 
is the particle position and k B T is the thermal energy, we 
arrive at the dimensionless Langevin equation describing 



2 



the overdamped motion of particles, i.e., 



dr 

d7 



F(r) +£(*), 



(2) 



with the Gaussian random force £ = £ y , £ 2 ) T obeying 
(&(*)) = 0, (&(*)£/(«)) = 25y£(*-s); i, j being either x, y 
or z. Although formally Eq. ^ is stated for a quiescent 
liquid, the case of driven solvent can be treated via the 
term V x *(r) entering F(r), see below. 

Let us next develop a generalized FJ theory for the 
formulated setup in Eq. ([2]). We use the corresponding 
Smoluchowski equation for the joint probability density 
function (PDF) P(r, t) |M1 



d t P(r, t) = - V [F(r)P(r, t)] + V 2 P(r, t) . 



(3) 



supplemented by the no-flux boundary conditions at the 
walls and the requirement of periodicity along the chan- 
nel (here x). In the spirit of the FJ approach, we 
perform the long-wave analysis [4T)lfi2"] in the dimen- 
sionlcss geometric parameter e = (Ail — Au>)/L <C 1 
[4151 21], which characterizes the deviation of the cor- 
rugated walls from flat structures, i.e., e = 0; Au and 
AQ denote the minimal and maximal channel widths, 
respectively. Upon re-scaling the transverse coordinate 
y sy, the profile functions become oj±(x) —> eh±(x), 
and *S? — > (e e \l/ 2 ) . Expanding as well the 

PDF in a series in even powers of e, we have P(r,t) — 
P (r, t) + e 2 P x {v, t) + 0(e 4 ), $(r) = $„(r) + 0(e 2 ), and 
similarly for ^(r), i — x,y,z. Substituting this ansatz 
into Eq. ^ and observing the boundary conditions, we 
obtain a hierarchic set. For the steady state, we find 
Po(r) = g(x,z) exp(— $ (r)), where g(x, z) is obtained 
in the order 0(e 2 ). 

With the conditions that (i) the x-component of 
V x 1 i r o( r ) is periodic in x with unit period and (ii) its z- 
component vanishes at z = 0, H, the stationary marginal 
PDF, P (x) = lim^ J^+g Ay J Q H dz P(r, t), yields 



P (x) =T- 1 I(x) 



(4) 



Here, I(x) = e"^) f* +1 dx'e^ x '\ 1 = dxl(x), and 
-P(x) is the generalized potential of mean force, reading 



F(x) = - In 



dy J dze-* ^ 

h-(x) 
h+{x') H 



(5) 



Jdx' J dy J dz (V x * ) x P cq (y, z\x') . 



h-(x') 



with P cq (y, z\x) = e-*»« / dy J Q H dz e -*°« . We 

reveal that .F(x) comprises the usual entropic contri- 
bution (the logarithmic term) [351 [M] caused by the 



non-holonomic constraint stemming from the boundaries 
[451 146j and the newly contribution, the part stemming 
from \I/o> which is associated with the conditional aver- 
age of the x-component of divergence-free forces exerted 
on the particle weighted by its equilibrium conditional 
PDF P cq (y,z\x). In the absence of *, Eqs. @ and ^ 
reduce to the commonly known result of the Fick-Jacobs 
approximation [24) [36], 

The kinetic equation for the time-dependent marginal 
PDF Po(x, t), with the steady-state solution in Eq. Q, 
is the generalized Fick-Jacobs equation, which reads 



Q , . d 
dt Po{x ' t)= dx 



dT(x) 
dx 



Po(x,t) 



+ J^oOM)- (6) 



We evaluate the stationary average particle current by 
use of well-known analytic expressions [36l [47] , to yield 



(x) = I~ 1 



(7) 



The effective diffu- 

v2\ 



wherein AT = F(x + 1) — T(x). 

sion coefficient D e ff = hmt^oo((^ 2 (^)) — ( x (t)) )/(2t) (in 
units of the bulk diffusivity, Dq = ksT / (Qir-qR)) is calcu- 
lated via the first two moments of the first passage time 
distribution, see Eq. (17) in Ref. [35], leading to 



D, 



off 



r-3 



d.r 



dx' e^-^')/ 2 (x). 



(8) 



In order to elucidate this result, we apply it next to 
Brownian motion under the influence of both, an exter- 
nal constant bias with magnitude / in x-direction, re- 
sulting in $(r) = —fx, and to the Stokes' drag force 
caused by the difference between the particle velocity r 
and the solvent flow field u(r) = V x \P(r). This implies 
a one-way coupling between the solvent and the parti- 
cles, when only the particle dynamics is influenced by 
the fluid flow but not vice versa |3Sj, as ensured by the 
adopted assumption of a dilute suspension. Accordingly, 
the particles dynamics is described by Eqs. ([!]) and 

Having mainly microfluidic applications in mind we 
shall focus on a slow pressure-driven steady flow of an 
incompressible solvent, determined by the dimensionless 
Stokes or "creeping flow" equations [12 [33] ; 

Vp(r) = V 2 u(r) , V-u(r) = 0, (9) 

being valid for small Reynolds number Re = 
pL 2 /(r]r) <C 1. Here, the flow velocity u= (u x ,u y ) T 
and the pressure p(r) are measured in the units of L/t 
and rj/r, respectively. We require that u obeys period- 
icity, u(x, y) = u(x + l,y), and the no-slip boundary 
conditions, u(r) = 0, Vr e channel wall. The pressure 
satisfies p(x + 1, y) = p[x, y) + Ap where Ap is the pres- 
sure drop along one unit cell. 

As the channel's height is much larger than all 
other length scales, we focus on the two-dimensional 
flow of incompressible fluid. Applying the curl to 
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both sides of first relation in Eq. Q eliminates p(r), 
yielding the biharmonic equation V ^(x, y) — for 
the stream function *&(x,y), = ^>(x,y)e z . Then, 
the components of the flow velocities are given by 
Ux = d y "i> and u y — —d x ^. With the above scaling, 
y — > ey, — > e^>, solving the biharmonic equa- 
tion = dy^o(x, y) + 0(e 2 ), and satisfying the no-slip 
boundary conditions, d y ^Q = at y = h±(x), and the 
conditions specifying the flow throughput, ^> — at 
y = h-(x) and * = - Ap j '(12("H~ 3 '(x)) x ) at y = h + (x) 
[4"5l [50] , we find in leading order the result 



Ap [y- h.(x)} 2 [3h + (x) - L(x) - 2y] 
12 ~ H*(x) (H-*(x)) x 



(10) 



where W(x) = h + (x) — h_(x) is the re-scaled local width 
and (•}.,,= Jq ■ dx denotes the average over one period of 
the channel. 

To elucidate the intriguing features caused by the 



divergence-free force based on Eq. ( 10 ) and its interplay 
with the constant bias we consider a reflection symmetric 
sinusoidally-shaped channel [4"31 144) . cf. Fig.JI] 



U!± (x) = ± 



An + Aw Afl - Aw 



cos (2nx) 



(11) 



Note that in the limit of the straight channel, 
S := Auj/An = 1, Eq. ^ yields the Poiseuille flow, 
u x = Ap [y 2 — (A£7/2) 2 ]/2 and u y = 0, between two 
plane walls at y — ±Afi/2. 

We next investigate the dependence of the transport 
quantities, such as the average particle velocity (x) and 
the effective diffusion coefficient D c g on the force mag- 
nitude / and the pressure drop Ap, which control the 
curl-free and the divergence-free contributions in Eq. (JlJ . 

Figure [2ja) depicts the impact of the pressure drop Ap 
on the mean particle velocity (£), cf. Eq. Q, for different 
values of external bias /. Only at / = 0, (x) = (x) 
point symmetric with respect to Ap, where 



is 



/Ap 



AAp(Aco) 2 VS 
'3(1 + 6) (3 + 25 + 35 2 ) 



f = 0. (12) 



The behavior changes drastically for / ^ 0. For 
Ap < with \Ap\ ^> 1, u x and / are both positive, the 
Stokes' drag dominates over the constant bias and thus 
(x) w (x) Ap oc — Ap. The increase in Ap results in a sys- 
tematic crossover from the flow-driven transport to bi- 
ased entropic transport. We observe a broad range of 
|Ap| with the width cx / in which the presence of the 



flow is insignificant, yielding (x) 
/ 3 + 4tt 2 / 



(x) 



f 



f 2 + 2ir 2 (V6 + l/V5) ' 



Ap = 0. 



(13) 



Note that for Ap > 0, the solvent flow drags the par- 
ticles into the direction opposite to the external force 
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Figure 2. (Color online) Comparison of Brownian dynam- 
ics simulations (markers) based on Eqs. |TJ and Q with the 
Fick-Jacobs approxi mati on (lines) for a corrugated channel 
with the profiles Eq. {n) for Afi = 0.5 and Au = 0.1. Panel 
(a): Mean particle velocity (x) versus pressure drop Ap for 
different force magnitudes /. The solid line corresponds to 



/ = 0, Eq. 1 12 1, the dashed lines represent Eq. |7| for / / 0. 
Panel (b): The depe nde nce of / cr on Ap cr for different S is 
depicted (lines: Eq. (14|). Panel (c): Effective diffusion D B g 
as a function of Ap for different / (lines: Eq. H) and the 
horizontal dash-dotted line corresponds to D e g = Do- Panel 
(d): Stationary joint PDF P(x,y) (color coding) obtained via 
direct simulation of Eqs. QTp , (|3J) , an d Q and force field F(r) 
(lines) for / = 10 2 and Ap = 6.5 ■ 10 4 , showing hydrodynam- 
ically enforced entropic trapping. 



(u x < and / > 0) and with increasing growth in Ap a 
sharp jump of (x) from positive to negative values occurs. 
Although strong nonvanishing local forces fe x + u(x, y) 
are acting on the particles, there exists a critical ratio 
( f/Ap) ci such that (x) = 0. As follows from Eq. Q, this 
occurs when F(x + 1) — F(x) = AF = 0, yielding for the 
critical ratio 



Ap 



1 (W(x)- 1 ) o 
12 (W(x)- 3 ) x 



2 An 2 s 2 



3 (3 + 2S + 3S 2 ) 



(14) 
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being solely determined by the channel geometry, see 
Fig. [2jb) . Here, W(x) denotes the local channel width, 



W(x) 



-(x) —oj-{x). Upon further increasing Ap, the 



flow-induced force starts to dominate over the static bias 
/ again and (A) s» (±) Ap cx - Ap. 

The role of Ap and / on the effective diffusion coeffi- 
cient D e ff is presented in Fig. [2^c). In the purely flow- 
driven case, / = (squares), D eS = 2y/S/(l + S) @3j for 
(A) < 1, i.e. small \Ap\. It exhibits so termed Taylor- 
Aris dispersion [ST] [52]; i.e., D cS cx (Af2 (±)) 2 /192 when 
(A) >■ 1, i.e. large \Ap\. 

In the limit of a resting fluid, Ap = 0, such that solely 
static bias induced transport occurs, the effective diffu- 
sion D e e exhibits the known bell shaped behavior as a 
function of / 34J. An intriguing effect emerges when 
the Stokes drag (u) and the external force (fe x ) exerted 
on the particle start to counteract, when u x cx — Ap 
and / are comparable, but of opposite signs. In this 
case, their superposition, F(r) = fe x + u, contains vor- 
tices and stagnation points, leading to hydrodynamically 
enforced entropic trapping (HEET). At a given / and 



Ap determined by Eq. ( 14 ) , yielding a vanishing particle 
current, D c g displays an abrupt decrease and is several 
orders of magnitudes smaller than the bulk value. Al- 
though the particles experience continuous thermal fluc- 
tuations, they exhibit long residence times in the domains 
of strong accumulation where the force field pushes the 
particles towards the channel wall, see Fig. [2jd) . This 
HEET-cffcct becomes more pronounced for larger / cr and 
Ap cr , resulting in a more localized particle distribution 
or, equivalently, larger depletion zones. This clarifies why 
the minimum of D e g decreases with the growth in /, see 
Fig. [2j leading to a stiffer trap. 

HEET offers a unique opportunity to efficiently sep- 
arate particles of the same size based on their different 
response to applied stimuli, e.g., to sift healthy cells from 
deceased and dead cells [HI [531 IM] • Even small distinc- 
tions in the response can be used to trap healthy cells 
and achieve opposite transport directions, cf. Fig. [3j for 
the deceased and dead cells by tuning / at a fixed Ap (or, 
equivalently, Ap at a fixed /) such that //Ap is close to 
the value given by Eq. ( fl4| ). We stress that the corru- 
gation of the channel, 8 ^ 1, is a crucial prerequisite for 
the function of an entropic sieve. For straight channels, 
5 = 1, the force field F(r) lacks vortices, the latter being 
responsible for particle accumulation. As a result, the 
effective diffusion coefficient is bounded from below by 
the value of bulk diffusivity and HEET fails. Thus, the 
Peclet number | (A) \/D g, which qualifies the transport 
of the objects, is strongly reduced compared to channels 
with finite corrugation, 5 ^ 1, cf. inset in Fig. [3] 

In conclusion, we generalized the Fick- Jacobs approx- 
imation for the most general force acting on the parti- 
cle, Eq. ([!]), which can involve both, the curl- free and 
the divergence-free components. Focussing on a typical 
corrugated channel geometry, we put forward an effec- 




Figure 3. (Color online) Snapshot of marginal PDF P(x, t) 
at t — 100 for different force strengths / in units of / cr = 100 
in a corrugated channel with Auu = 0.1. The width of P(x,t) 
is several magnitudes smaller compared to the case of un- 
bounded geometry cx Do t. Inset: Peclet number | (x) \/D e s 
versus ///cr for 8 = 0.2 and 5 = 1 (straight channel). The 
maximum width Afl — 0.5 and Ap = 6.5 • 10 4 are kept fixed. 



tive one-dimensional description involving the potential 
of mean force, which along with the commonly known 
entropic contribution in presence of a constant bias, ac- 
quires a qualitatively novel contribution associated with 
the divergence-free force. 

The analysis of particle transport caused by the coun- 
teraction of a pressure-driven flow (presenting the case of 
a divergence-free force) and a constant bias of strength 
/ pointing in the opposite direction, ensues the intrigu- 
ing finding that the mean particle current can identi- 
cally vanish despite the presence of locally strong forces. 
Being accompanied by a significant suppression of diffu- 
sion, thus being robust against thermal fluctuations, this 
purely entropic effect of strong particle accumulation, in- 
duced by the corrugation of the channel, yields a selective 
hydrodynamically enforced entropic trapping, which can 
be utilized to separate particles of same size. The the- 
oretical predictions here are in excellent agreement with 
the results obtained from numeric simulations. Note that 
our methodology admits the situation of a driven solvent; 
alternatively, similar effects can be expected in a resting 
solvent with nonvanishing divergence-free forces. 
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